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ANOMALOUS  DISPERSION  IN  GASES  DERIVED  FROM  THE  OPTICAL  DEPTH 
Theoretical  treatment;  line  by  line  calculations 


1  INTRODUCTION 

This  report  describes  the  high  resolution  wave  numoer  dependence  of 
the  refractive  index  in  gases.  The  wave  number  dependence  of  the 
refractive  index  is  caused  by  the  ability  of  a  gas  to  absorb  and  emit 
electromagnetic  radiation.  Knowledge  of  the  spectral  absorption  (or 
the  optical  depth)  is  the  basis  for  refractive  index  calculations  in 
this  report. 

The  spectral  absorption  can  be  calculated  from  the  quantity  optical 
depth.  The  spectral  behavior  of  the  optical  depth  is  mainly  described 
by  the  spectral  line  shape  function.  The  correct  line  shape  function 
is  not  known,  but  the  approximations  of  Voigt,  Lorentz  and  Doppler  are 
normally  used  for  calculation  of  spectral  absorption.  These  approxi¬ 
mations  will  also  be  used  in  this  report  for  the  refractive  index 
calculations1.  They  give  a  correct  description  of  the  line  shape  near 
the  resonance  wave  number  but  far  away  from  resonance  they  are  not 
correct.  This  effect  has  greater  implications  on  the  accuracy  of  cal¬ 
culated  refractive  index  than  calculated  absorption  due  to  differences 
in  the  spectal  behavior. 

Line  by  line  calculation  of  absorption  at  a  certain  wave  number  is 
normally  performed  by  only  taking  into  account  the  influence  from 
transitions  "nearby"  the  given  wave  number.  The  contribution  to  the 
absorption  from  the  "wings"  of  distant  transitions  is  neglected.  This 
approximate  calculation  of  absorption  gives  fairly  good  results 
because  the  line  shape  function  falls  off  rapidly  (~(o-Oj)-2)  with 
the  distance  from  the  resonance  wave  number  (o^). 


*A  correct  line  shape  function  would,  among  other  effects,  have  to 
take  the  effect  of  line  coupling  into  account.  This  effect  will  not 
be  discussed  in  this  report. 
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Line  by  line  calculation  of  refractive  index  leads  to  a  somewhat  dif¬ 
ferent  situation  with  respect  to  accuracy  than  encountered  in  the  cal¬ 
culation  of  absorption.  The  reason  is  that  the  contribution  to  the 
refractive  index  from  the  "wings"  of  distant  transitions  is  much 
greater  due  to  the  slower  fall  off  (~(o-o-j)_l)  with  the  distance  from 
the  resonance  wave  number.  This  means  that  calculation  of  the  refrac¬ 
tive  index  at  a  given  wave  number  requires  that  the  contribution  from 
the  "wings"  of  a]_l_  resonance  transitions  are  taken  into  account.  When 
dealing  with  atmospheric  refractive  index  especially  the  contribution 
from  the  very  strong  UV-transitions  has  to  be  taken  into  account  both 
in  the  visible,  infrared  and  y-wave  part  of  the  spectrum. 

A  computation  of  the  refractive  index  as  outlined  above  is  not  treated 
in  this  report  although  the  theoretical  results  of  Chapter  4  are 
suitable  for  such  computation  provided  the  Voigt  or  Lorentz  line  shape 
function  is  used.  The  result  from  such  computation  will  have  limited 
value  due  to  the  validity  of  the  distant  "wings"  of  the  line  form 
functions.  The  computations  in  this  report  is  performed  by  only  tak¬ 
ing  into  account  the  contribution  to  the  refractive  index  from  reso¬ 
nance  transitions  in  a  limited  wave  number  interval  centered  around 
the  wave  number  of  interest  (i  e  a  procedure  analogous  to  the  one 
described  for  calculation  of  absorption).  The  result  from  these  com¬ 
putations  will  therefore  only  give  the  local  wave  number  dependence  of 
the  refractive  index;  more  commonly  named  the  anomalous  dispersion  due 
to  resonance  transitions.  This  quantity  is  in  this  report  called  the 
dispersion  and  is  denoted  An. 

For  the  atmospheric  refractive  index  an  empirical  formula  has  been 
developed  by  Edlen  (1)  based  on  measurements  in  the  visible  part  of 
the  spectrum.  This  formula  with  corrections  for  atmospheric 
temperature,  humidity  and  pressure  is  also  widely  used  for  calculation 
of  atmospheric  refractive  index  in  the  infrared,  y-wave  and  even 
radio-wave  part  of  the  spectrum.  The  accuracy  of  these  results  is 
however  questionable  since  the  basic  measurements  are  only  performed 
in  the  visible.  It  is  however  reasonable  to  believe  that  Edlen’s  for 
mula  expresses  the  total  contribution  to  the  refractive  index  from  the 
"wings"  of  distant  strong  transitions  especially  the  very  strong 
UV-transitions  In  the  atmospheric  spectrum.  The  coarse  (low 
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resolution)  wave  number  dependence  of  atmospheric  refractive  index  is 
then  described  by  Edlen's  formula  whilst  the  fine  structure  (high 
resolution)  wave  number  dependence  is  described  by  the  dispersion 
derived  in  this  report.  A  sum  of  these  two  quantities  is  then  an 
estimate  of  the  total  atmospheric  refractive  index.  The  contribution 
fn  this  sum  from  the  dispersion  is  small  in  the  visible  and  infrared 
part  of  the  spectrum.  In  the  wave  number  region  below  400  cm"1  the 
dispersion  is  important  near  the  resonance  transitions.  This  situation 
will  be  demonstrated  and  discussed  in  Chapter  6. 

The  calculations  in  this  report  is  motivated  by  the  fact  that  there 
already  exists  a  computer  code  performing  line  by  line  calculations 
of  the  optical  depth.  The  code  is  named  Fast  Atmospheric  Signature 
Code  (FASCOD),  and  is  developed  by  the  U  S  Air  Force  Geophysics 
Laboratory  (AFGL) .  In  Chapter  5  it  will  be  shown  how  line  by  line  cal 
culations  of  the  dispersion  can  be  performed  within  the  framework  of  a 
modified  FASCOD. 

This  report  consists  of  two  main  parts:  A  theoretical  development  of 
the  expression  for  the  dispersion  (Chapter  2-4)  and  an  implementation 
of  the  theoretical  result  in  FASCOD  (Chapter  5-6). 

In  Chapter  2  it  is  shown  how  the  dispersion  can  be  expressed  in  terms 
of  the  optical  depth  by  means  of  the  linear  susceptibility  and  the 
Kramers-Kronig  relations.  The  formulation  of  the  optical  depth  used 
in  this  paper  is  discussed  in  Chapter  3.  An  explicit  expression  for 
the  dispersion  is  calculated  in  Chapter  4. 

A  theoretical  development  analogous  to  the  one  described  in  this  report 
is  given  by  Tomiyama  (2).  The  theoretical  result  in  this  report  has 
however  a  more  general  validity  since  no  approximations  are  needed 
during  the  development.  This  gives  an  exact  expression  for  the  dis¬ 
persion  when  the  Voigt,  Lorentz  or  Doppler  line  shape  functions  are 
used. 

In  Chapter  5  it  is  shown  how  the  expression  for  the  dispersion  can  be 
simplified  to  facilitate  implementation  within  the  framework  of 
FASCOD. 
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Results  from  calculated  amtospheric  dispersion  in  various  wave  number 
regions  are  presented  and  discussed  in  Chapter  6.  The  calculated 
results  are  compared  with  the  measurements  of  Liebe  (3). 


2  DISPERSION  AND  OPTICAL  DEPTH 

This  chapter  describes  how  the  dispersion  can  be  expressed  in  terms  u< 
the  optical  depth  by  means  of  the  linear  susceptibility  and  the 
Kramer s-Kronig  relations. 

The  linear  susceptibility,  X(u>),  is  a  complex  quantity  describing  the 
relationship  between  an  external  electric  field,  ?,  and  the 
polarization,  that  is  induced  by  the  electric  field.  When  the 
medium  is  isotropic  the  linear  susceptibility  is  a  scalar,  and  the 
polar izat ion  is  given  by 

f5  =  x(w)f  =  (x'(w)  +  lX”(«))t  (2.1) 

where  x'(w)  is  the  real-  and  x"(w)  the  imaginary-part  of  the  suscepti 
bility.  u  is  the  angular  frequency. 

In  this  context  we  shall  assume  a  plane  wave  propagation  for  the 
polarization  ? 

?  .  »  f  ei (wt-KZ)  ,  c>c  (2.?) 

p  is  the  complex  magnitude  and  k  =  x'+ix"  is  the  complex  propagation 
constant  with  real  part  k'  and  imagninary  part  x". 

The  large  scale  macroscopic  phenomena  of  absorption,  emission  and  dis¬ 
persion  of  light  are  the  results  of  transitions  between  the  various 
atomic  or  molecular  energy  levels.  The  equation  of  motion  for  the 
electric  and  polarization  fields  in  an  isotropic  medium  are  derived  in 
Pantell  (4).  When  a  plane  wave  solution  of  the  form  (2.2)  is  assumed, 
the  equations  of  motion  gives  the  following  expression  for  the  propa 
gation  constant,  x 
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'2.3) 


In  Equation  (2.3)  n  is  the  contribution  to  the  propagation  constant 
from  transitions  outside  the  frequency  interval  considered  for  deter 
mination  of  the  susceptibility  (i  e  the  transitions  that  are  not  taken 
into  account  in  the  calculations).  The  value  of  n  is  equal  to  unity 
if  the  contribution  from  aJJ_  transitions  in  a  gas  is  accounted  for  in 
the  the  determination  of  the  susceptibility. 


Equation  (2.3)  is  soivea  by  taking  the  square  root  and  assum¬ 
ing  (X(w)/n2 )  <<  1  so  that  the  expansion  [i+x(w)/nJ]^  =  l+X(w)/2r,2 
ho'ds  (which  is  the  case  for  other  than  very  strong  transitions).  I  he 
result  for  the  complex  propagation  constant  is 


x  -  -  '  +  i x "  =  [ ] 

The  real  and  imaginary 
by ; 


part  of 


2qc 


the  propagation  constant  is  then  given 


(2.5) 


and 


K" 


u>X"(Q) 

?nc 


(2.6) 


Since  the  optical  depth,  k(w),  is  defined  with  reference  to  the 
absorption  of  radiant  power  it  can  be  shown  (4)  that  the  imaginary 
part  of  the  propagation  constant  is  related  to  the  optical  depth  by 


k(u)  =  -  2k"  = 

v  nc 


(2.7) 


The  refractive  index,  n(w),  in  the  medium  is  defined  as  the  ratio 
between  the  phase  velocity  in  empty  space,  c,  and  the  phase  velocity 
in  the  medium,  vf.  By  definition  the  phase  velocity  in  the  medium  is 


(2.8) 


The  relation  between  the  refractive  index  and  the  real  part  of  the 
susceptibil ity  is  then 
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n(u) 


The  dispersijn  as  defined  in  this  report  is 


(2.9) 


&n(w)  = 


(2.10) 


i he  Equations  (2.7)  and  (2.10)  connecting  the  imaginary  and  real  part 
of  the  susceptibility  to  respectively  optical  depth  and  dispersion, 
are  the  fundamental  equations  in  this  report. 


By  us’ng  the  connections  between  angular  frequency  and  wave  number 
(u>  2ttco)  these  equations  can  be  rewritten 


,Vi(o) 


L 

2n 


X'  (o) 


(2.H) 


Mo)  --  -  —  xN(0)  (2.1/; 

the  purpose  of  this  chapter  is  to  express  the  dispersion  in  Equation 
(2.J1)  by  the  optical  depth  in  Equation  (2.12). 

As  stated  by  the  Kramer-Kronig  relations  (4)  the  real  and  imaginary 
part  of  the  susceptibility  are  Hilbert  transformed  pairs.  Mathemati 
caily  this  can  be  expressed 


x *  (0) 

=  -  h[  x"  ( o )  ] 

(2.13) 

X"(o) 

=  Hi  X* (o) J 

(2M4) 

The  Hilbert  transform,  H( f ( x )  ] ,  of  a  function  f(x)  is  defined  by 


Hlf(x)]  =  l  PP  /  dx*  (2.  IB) 

«e» 

where  PP  means  that  the  principal  part  of  the  integral  is  to  be 
evaluated. 

From  Equations  (2.12)  and  (2.13)  the  relationship  between  the  real 
part  of  the  susceptibility  and  the  optical  depth  is 


X  '  ( o ) 


(2.16) 


iO 

=  -  "I-  JT5  1(011  ’  5i  "1^1 

"he  Equations  (2.11)  and  (2.16)  then  establish  the  following  relation 
between  the  dispersion  and  the  optical  depth 

An(o)  =  ^  j 

The  relationship  expressed  by  Equation  (2.i7)  is  the  basis  for  the 
description  throughout  the  rest  of  this  report. 

in  order  to  evaluate  the  Hilbert  transform  in  Equation  (2.17)  the 
explicit  wave  number  dependence  of  the  optical  depth  must  be  knovn. 
This  p^oMem  will  be  addressed  ’n  the  following  chapter. 


I  THE  OPTICAL  DEPTH  FORMULATION 

Ttiis  chanter  contains  a  formulation  of  the  optical  depth  which  takes 
into  account  the  molecular  processes  of  absorption  and  stimulated 
emission.  The  explicit  wave  number  dependence  oi  the  optical  depth  ^ 
also  described  both  for  the  Voigt.  Lorentz  and  Doppler  line  shape 
functions.  The  aim  o *  this  chapter  is  to  define  that  part  of  the 
optical  depth  which  has  to  be  Hilbert  tranfo'med  in  order  to  calculate 
the  dispersion  and  to  define  the  Voigt,  Lorentz  and  Doppler  dispersion 
functions. 

The  optical  depth  formulation  of  Van  Vieck  and  Huber  (5),  will  be  used 
in  this  paper.  In  this  .ormulation  the  optical  depth  k(ol  of  a 
mixture  of  different  gases  can  be  expressed  in  the  following  form  (61 

o  tanh  ,'|2)  l  W(m . )  5,(1)  [f(o,oJ  +  f(-o,o.)]  (3.1) 

i  j  J  1  1 

-  Number  density  of  molecular  species  mj  (molecule/cm5) 

-  Line  strength  of  *he  1‘th  transition  at  the  temperature  1 
(cm2/molecule) 


k  ( o )  = 

where 

W(TTIj) 

Si(T) 


I 


f(0,0  ) 
0 
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Line  shape  function  (cm) 

Wave  number  (cm-1) 

Resonance  wave  number  for  the  i '  th  transition 

=  hc/kT;  where  T  ic  temperature  and  hc/k  is  the  second 
radiation  constant  (1.4388  cm*k) 


The  Tine  strength.  S-j(T),  in  Equation  (3.1)  has  a  specific  value  for 
ei  h  resonance  transition  belonging  to  the  different  molecular  species 
in  the  gas  mixture.  The  temperature  dependence  of  the  line  strength 
is  given  by  the  assumption  of  thermal  equilibrium;  i  e  a  Maxwell- 
8nlt.zmann  .istribution  of  molecules  in  the  excited  states  and  the 
associated  partition  functions. 

This  formulation  of  the  optical  depth  satisfies  three  important  ^rop 
erties  independent  of  the  iine  shape  function.  These  properties  are 
satisfied  even  if  the  fluctuation  dissipation  theorem  for  detailed 
radiation  balance  is  not  (5).  The  properties  satisfied  by  Equation 
(3,1)  are: 

a)  i he  generalized  Nyquist  theorem;  i  e  transition  strengths  arc 
preserved,  and  radiation  balance  between  emission  and  absorption  is 
satisfied  in  thermal  equilibrium  (7). 

b)  The  kramers-Kronig  relations,  (5)  and  (8). 

c)  The  f  sum  rule;  which  's  equivalent  to  the  validity  of  the  familiar 
formula  (qp-pq)  =  lh/2i?  nf  matrix  quantum  algebra  (5). 

This  means  that  Equation  '3.1)  should  be  valid  for  aH  wave  numbers 
provided  the  dipole  approximation  is  valid  and  a  classical  despribticn 
of  the  radiation  field  is  adequate.  The  accuracy  of  the  caloulated 
optical  depth  is  only  dependent  of  the  quality  of  the  approximate  line- 
shape  function. 

In  Chapter  2  the  connection  between  dispersion  ana  optical  depth  is 
expressed  by  Eauai.ion  (2.17).  By  introducing  the  optical  depth  in 
Equation  (3.1)  the  dispersion  can  be  expressed  ty 


-  k  I  Wfi.i.  )S.(r)[F(a.a1) 
ij  J 


-*  G(o,o.)] 


An(o) 


(3.2) 
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■where 


f  ( 0 , 0  1  ) 

-  h[  t.anh 

(|?i 

•  f(o,o1)  ) 

(3.3) 

G  ( o .  c  , ) 

i 

-  h[ tanh 

(|2) 

’  f(-o,oi)] 

(3.4) 

H [  j  is  the'Hil&ert  tranform. 


To  calculate  the  Hilbert  tranforms  described  by  Equations  (3.3)  and 
(3.4)  it  is  necessary  to  introduce  an  explicit  expression  for  the  line 
shape  function,  f(o,Oi). 

i  r>  order  to  obtain  a  theoretical  result  for  the  dispersion  which  con- 
tain  the  influence  of  both  homogeneous  and  inhomogeneous  line  broaden¬ 
ing  mechanisms  the  Voigt  line  shape  function  is  used  for  the  theoreti 
cai  calculations  in  this  report.  The  Voigt  line  shape  is  a  convolu¬ 
tion  of  the  Lorent2  line  shape  describing  homogeneous  broadening  due 
to  molecular  collisions  and  the  Doppler  line  shape  describing 
inhomogeneous  broadening  due  to  the  distribution  of  molecular 
velocities.  The  Voigt  line  shape  function  is  described  as 


V0-0,)  -  *  ~  i  l  ro-~ exp  [-In?  (-—V  ]dx  (3.5) 

o  -  ®  a 

The  Lorentz  and  Doppler  line  shape  functions  are  defined  by 


,(o.oi ) 


l  a 

TI  (O-O^5  +  <P 


(3.6) 


1  0-0.  * 

—  exp  [  -  l  n2  (  --■--)  |dx 
d  a 


where 

a  Lorentz  half  width  (HWHM) 

ad  -  Doppler  half  widt  (HWHM) 


(3.7) 


The  half  widths  are  dependent  of  molecular  species  and  the  thermo¬ 
dynamic  situation  in  the  gas  mixture.  The  Doppler  line  width  is  also 
proportional  to  the  resonance  wave  number.  At  standard  atmospheric 


temperature  and  pressure  the  Lorent2  line  width  is  normally  much 
greater  than  the  Doppler  line  width. 


1  he  Voigt  line  shape  in  Equation  (3.5)  is  easily  obtained  by  a  convo¬ 
lution  of  the  Lorentz  and  Doppler  line  shapes  in  Equations  (3.6)  and 
(3.7)  (regarding  the  wave  number  difference  o-o-j  as  free  variable) 
followed  by  a  simple  change  of  variable. 

V0’0,)  *  *  fd(°-oi) 


f® 

?  /  f2(o-o1  y)  •  fd(y)dy 


--  f  fa(o-x)  •  f(J(x-o1)dx 


(3.8) 


By  using  the  notation  in  Equations  (3.6)  and  (3.7)  the  last  expres¬ 
sions  in  Equation  (3,8)  can  formally  be  written 


+® 

*  J  f4(°’x)  •  fd(x,o.)dx 


(3.9) 


The  notation  in  Equation  (3.9)  means  that  the  variable  indicated  (with 
x)  in  the  Lcrentz  and  Doppler  line  shape  functions  shall  be  substi¬ 
tuted  by  the  integration  variable  prior  to  integrations.  This  type  of 
notation  is  convenient  at  a  later  stage  in  this  report. 

The  Voigt  line  shape  functions  is  within  the  approximations  of  Lorent? 
and  Doppler  line  broadening  mechanisms  a  reasonable  description  of  the 
line  shape  (to  a  certain  distance  from  the  line  center)  for  most  values 
of  temperature  and  pressure  in  a  gas  mixture.  However  at  two  extreme 
situations  the  Voigt  formula  can  be  simplified  since  (9) 

11m  fd(o,o.)  =  5(0-0.)  (3.10) 

V° 

lim  f  (o,o.)  =  S(o-oJ  (3.11) 

a->0  *  1  1 

where  5(*)  is  the  Dirac  5-function. 
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The  Equations 
for  the  Voigt 


(3.10),  (3.11)  and  (3.9)  leads  to  the  following  limits 
line  shape  function  when  the  half  widths  tends  to  zero 


Mm  f  (o.o^  = 

a,~*0 

a 


Mm  f  (o,o.)  = 

a->0 


(3.1?) 

(3.13) 


Equation  (3-12)  describes  a  situation  where  homogeneous  (Lorentz) 
broadening  completely  dominates  over  inhomogeneous  (Doppler) 
broadening,  i  e  the  high  pressure  regime.  Equation  (3.13)  describes 
the  low  pressure  regime  when  >>  a.  Since  the  natural  line 
broadening  due  to  spontaneous  processes  is  a  homogeneous  broadening 
mechanism  Equation  (3.13)  is  only  a  good  approximation  for  the  low 
pressure  regime  when  the  Doppler  line  width  is  much  greater  than  the 
natural  line  width;  >>  an.  The  natural  line  width  is  given  by 


<■„  -  (cinM 

where  xn  is  the  natural  lifetime. 


(3.14) 


To  make  this  theoretical  development  valid  for  both  the  high  and  low 
pressure  regime  the  Voigt  line  shape  function  is  used  in  Equation 
(3.2).  The  fact  that  fv(-o><Ji)  -  fv(°  -o-j )  (which  is  obvious  by 
change  of  variable  in  Equation  (3.5))  gives  the  following  expression 
for  the  dispersion 


An(o)  =  ^  X  W(m  )S  (T)[F  (o,0i)  4  F  (o,-o  )]  (3.15) 

i  j  j 

where 

Fy(°»o1 )  =  H[tanh  (|2)  fv(o.o.) J  (3.16) 

Fv(o,-0i)  =  H[tanh  (|2)  fy(o,-oi)]  (3.17) 


From  the  definition  of  the  Hilbert  transform  in  Equation  (2.15)  and 


15 


the  Voigt  line  shape  function  in  Equation  (3.5)  an  explicit  expression 
for  Equation  (3.16)  is 


-S'-T-r.  ■H|tanh(r>'l  i  rdrfwe><p  l-'nZ  'tt1’1'”" 


(3.18) 


]Jr  ‘  5'  _■[  -H'tanh(|2)-  £  ^ ]*exp  M"2  (~J)  ]d 


x-o. 


d 

(3,19) 


The  Hilbert  transform  in  Equation  (3.19)  is  given  by: 


F  2 ( o , x )  5  H| tanh  (|2)  • 


-  I  pp  f 

T1 


±  a  1 

n  (o~x) }+05  J 

tdnh  (f2-)^  do 


i pp  i 


i 


o'  -o 


tanh  (1^-) 


fa(o’ ,x)  do' 


(3.20) 


The  integral  in  Equation  (3.20)  will  be  calculated  in  Chapter  4  by 
analytic  continuation  and  complex  integration.  The  high  and  low  pres 
sure  limits  for  the  dispersion  are  found  from  the  limits  of  Equation 
(3.16)  using  (3.12)  and  (3.13) 


Pa(°.V  5  11m  Fv(0,0i)  =  l  PP*J  ^hro  tanh(^')  ♦  f^o'.o,)  do' 

V°  (3.21) 

+«* 

F.(o,o.)  =  lim  F  (o,a.)  -  ^  PP  /  tanh(^f’)  •  fd(o',o.)  do' 
a  1  a  -»0  v  1  ^  0  0  *  a  1  (3,22) 

The  functions  Fv(o,o1) .F^o.o^)  and  F^o.o^)  will  in  the  rest  of  this 
report  be  called  the  Voigt-,  Lorentz-  and  Doppler  dispersion  functions. 

A  comparison  of  Equations  (3.21)  and  (3.20)  shows  that  the  high  pres 
sure  limit  of  the  dispersion  Is  readily  calculated  line  by  line  from 
Equation  (3.15)  as  soon  as  the  integration  In  (3.20)  is  performed  and 
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the  variable  x  is  substituted  by  the  resonance  wave  numbers  o..  At 
intermediate  pressures  the  integration  described  by  Equation  (3.19) 
must  be  performed  before  a  line  by  line  calculation  of  the  dispersion 
can  be  evaluated. 

A  calculation  of  the  integral  in  Equation  (3.19)  is  difficult  within 
the  framework  of  FASCOD  therefore  the  [numerical  examples  given  in  this 
report,  will  mainly  concentrate  on  the  high  pressure  regime. 


4  THE  HILBERT  TRANSFORM  OF  THE  OPTICAL  DEPTH 

The  purpose  of  this  chapter  is  to  calculate  an  explicit  expression  for 
the  dispersion  based  on  the  optical  depth  formulation  in  Chapter  3  and 
the  Voigt  line  shape  function. 

The  development  in  Chapter  3  showed  that  the  following  Hilbert  trans¬ 
form  has  to  be  calculated  in  order  to  determine  the  dispersion. 


o.x) 


tanh  (^|~)  • 


I  a 

n  (o' -x) J+a2 


do ' 


(4.1) 


By  using  the  notation  F(o')  for  the  intergrand  in  Equation  (4.1)  the 
evaluation  of  the  principal  part  of  the  integral  is  by  definition 


,  o-e  +- 

F.(o,x)  «  i  11m  {  J  F(o' ) do  *  +  /  F ( o ' )do‘ }  (4.2) 

e-*o  -•  o+e 

To  be  able  to  evaluate  the  integral  of  Equation  (4.1)  a  continuation 
into  the  complex  plane  will  be  made.  By  defining  a  closed  contour  of 
integration  it  is  possible  to  use  the  residue  theorem  from  complex 
analysis  to  evaluate  the  integral.  The  contour  of  integration  which 
is  used  here  is  shown  in  Figure  4.1.  It  Is  called  C,  and  it  consists 
of  the  curves  r,f'  and  I. 
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k  Im  {O  } 
(J)  (2_n  _+  1)  n 

:  '  p" 


Figure  4,1  The  integration  contour  C  in  the  complex  plane  which  is 
used  in  the  evaluation  of  the  integral  in  Equation  (471) 
The  poles  of  the  integrand  of  Equation  (4.1)  are  indi- 
cated  by  the  symbol  O 

The  curve  r  in  Figure  4„1  is  a  semicircle  in  the  upper  halfplane  with 
radius  R.  It  meets  with  the  real  axes  at  the  points  o'  =-R  and 
o'  =  R,  and  it  intersects  the  imaginary  axes  at  the  point  o'  =  iR.  In 
appendix  A  it  will  be  shown  that  the  integrand  of  Equation  (4.1)  has 
infinitely  many  poles  on  the  imaginary  axes.  There  is,  however,  a 
finite  distance  unequal  zero  between  these  poles.  This  means  that  the 
curve  T  can  always  be  made  to  intersect  between  two  such  poles. 2  in 
the  limit  where  R-*~  the  integral  along  the  curve  T  will  tend  to  zero 
since  the  integrand  of  Equation  (4.2)  goes  to  zero  as  R~3  for  big  |  o’) 

The  curve  f  is  a  small  semicircle  with  radius  e  centred  in  o'=o. 

The  integral  along  this  curve  will  be  evaluated  for  the  limit  e-»0  in 
Appendix  A. 

The  curve  I  consists  of  two  Intervals  on  the  real  axis,  namely 
-R^o'$o-e  and  o+eso'<R,  where  e>0.  In  the  limit  R-*«  and  e->0  the  inte- 


2This  may  of  course  slightly  distort  the  value  of  the  Intergrand  on 
the  semicircle  r,  but  in  the  limit  R-»-  such  a  distortion  will  not 
change  the  value  of  the  integral  along  f. 
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gral  along  1  is  the  principal  part  of  the  integral  in  Equation  (4,1). 
By  means  of  the  residue  theorem  the  closed  contour  integral  is  given 
by 


J  F(o')d o'  =  J  F(or)do'  +  j  F(o')do’  +  j  F(o')do' 

c  i  r  r 

(4.3) 


„  .  [The  sum  of  the  residues  of 
^  [_ F ( a ’ )  inside  the  contour  C 

In  the  limit  R-+-  and  e-»o  the  result  is 


pp  J  F(o* )do' 


?iri 


The  sum  of  the  residues  of 
F(o')  inside  the  contour  C 


i im  J  F(o' )do 
e-*o  P 


(4.4) 

A  detailed  calculation  of  the  right  hand  side  of  Equation  (4.4)  is 
described  in  Appendix  Av  and  the  result  is  given  in  Equation  (A. 18). 


The  result  for  the  Hilbert  transform  described  by  Equation  (4.1)  is 


F 


(o,x) 


i 

i 


PP  /  F  (o  ’  )do  ’  -- 


-  x-o 
it  (x-o) J+aJ 


sinh(Bx) 

cosh(Bx) fcos (Ba) 


+  l  Im  {vP(^+i|i(x-ia))-T|;(Vi+i|i(x+ia))} 


1  a 
tt  (x-o)  !+aJ 


sin(Ba) 

cosh(Bx)+cos(Ba) 


{'H^+i|S(x-ia))+T|/(Vi+i|i(x+ia))! 


-  £  Re  {♦(VS+il^o)} 

The  function  t))(z)  is  the  digamma  function  defined  by 

»<*>  -  3z  C«rW)  . 

where  f(z)  is  the  gamma  function  (9). 


(4.5) 


(4.6) 


The  Voigt  dispersion  function  Fv(o,Oi)  defined  in  Equation  (3.19)  can 
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now  be  rewritten  by  Equation  (4.5) 


+® 


F.(0,0,)  ■  /'-¥  ■  ~  I  Ft(0.x) 

a  -°® 


x*o.  2 

exp  [  - 1  n2  (— — — i)  ]dx 
d 


(4.7) 


A  notation  equivalent  to  the  one  used  in  Equation  (3.9)  gives  the  foi 
lowing  compact  formula  for  the  Voigt  dispersion  function 


+• 

F  (0,0^  =  J  F  (o.x)  •  fd(x,o.)dx  (4.8) 

-  ao 

where  f d ( * , - )  is  the  Doppler  line  shape  function. 

The  high  pressure  limit  for  the  Voigt  dispersion  function  is  found 
from  Equation  (4.8)  and  (3.10) 

nm  F  (o,o  )  =  F  (o,o  )  (4.9) 

V° 

The  Lorentz  dispersion  function,  F4(o,Oj),  is  given  by  Equation  (4.5) 
when  the  variable  x  is  substituted  with  the  resonance  wave  number,  o,. 

The  low  pressure  limit  for  the  Voigt  dispersion  function  is  also  found 
from  Equation  (4.8)  when  the  limit  of  Equation  (4.5),  as  the  Lorentz 
half  width,  a,  tends  to  zero,  is  determined.  The  limit  of  Equation 
(4.5)  is 


1im  Fi(0’x)  =  i  zh  '  tanh(r) 

a-*o 

(4, JO) 

+  \  S(x-o)  •  Refy(f  +  if^x)  -  i|)(|  +  i|jjo)} 

The  low  pressure  limit  for  the  Voigt  dispersion  function  is  then 

lim  F  (o,o.)  =  ~  /  tanh  (§*)  fd(x,o,)dx  (4.11) 

a-*o  -«• 


Since  the  integrand  in  (4.11)  diverges  where  x=o  the  integration  must 
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be  performed  by  calculation  of  the  principal  part.  In  this  respect 
Equation  (4.11)  is  identical  to  Equation  (3.22)  defining  the  Doppler 
dispersion  function,  v.i.z 

Mm  F  (o,o  )  =  F  (o,oi)  (4,12) 

a-^o 

An  explicit  expression  for  the  Doppler  dispersion  function  will  not  be 
caicuaited  in  this  report. 

At  this  stage  it  is  useful  to  remind  oneself  that  the  dispersion, 
An(o),  given  by  Equation  (3.15)  has  the  following  form 

An(o)  --  V  W(m .)S  (T)[F  (o,o.)  +  F  (o.-o  )]  (4.13) 

i  j  J 

By  using  Equations  (4.5)  and  (4.9)  it  is  now  possible  to  give  an 
explicit  expression  for  the  dispersion  in  the  high  pressure  limit 


o.  -  o 


MO)  .  i-„  IW(V  3,(7)  i  •  7,(0,)^  Ti-T5fw  ■  T2(0.», 


1  V° 


+  11  (o^+o)  J-H3J  '  Tl^°i^  +  IT  ( 0  ^  +0  )  *  -KJ 2  "  T2^0,0i^ 


(4.14) 


The  functions  T^(*)  and  T2(*,*)  are  easily  identified  from  Equation 
(4.5).  To  establish  Equation  (4.14)  the  following  relations  are  used 


Ti(-Oi)  =  -T !  ( o  -j ) 

T2(o,-Oi)  =  T2(o,Oi)  (4.15) 

The  well  known  anomalous  dispersion  form  due  to  a  resonance  transition 
is  recognized  from  the  first  term  in  Equation  (4.14). 

To  Illustrate  the  wave  number  dependence  of  the  different  terms  in 
Equation  (4.14)  the  four  terms  (exclusive  the  Tj  and  T2  functions)  are 
plotted  separately  for  two  different  resonance  wave  numbers 
(0}  =  1  cm-1  and  o^  =  100  cm-1)  in  the  Figures  4.2  and  4.3.  The 
Lorentz  half  width  and  the  temperature  are  0.1  cm-1  and  296  K.  The 
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value  of  the  functions  Ti(oi)  and  T2 (o,Oi)0=0i  is  given  (on  the 
figures)  for  each  resonance  wave  number. 

Figure  4.2  illustrates  the  situation  when  the  resonance  wave  number  is 
1  cm*1  The  numbers  refering  to  the  different  curves  indicate  the 
term  number  in  Equation  (4.14). 


Figure  4.2  Wave  number  dependence  of  the  different  terms  in  Equation 
^4. 14)  for  resonance  wave  number  Qj  =  1  cm-1 
The  functions  Tj(a-j)  and  T2(oto-j  )0=0.  have  the  values 
indicated 


For  the  discussion  of  the  results  in  Figure  4.2  the  following  symbolic: 
representation  of  the  wave  number  dependent  part  of  Equation  (4.14)  is 
useful  (Assuming  T2(o,o-j)  is  a  slowly  varying  function  when  0=0 -j ) 

H  •  ®  *  T1(o1)  +  ©  •  T2(o.,0i)  +  ©  •  Tj (0. )  +  ©  *  T2(o . ,0 . ) 

(4.16) 

The  results  In  Figure  4.2  shows  that  the  peak  value  of  the  first  term 
in  Equation  (4.16)  is  dominating.  The  peak  values  of  the  second  and 
third  term  are  about  one  and  two  orders  of  magnitude  less  than  the 
first  term.  The  value  of  the  fourth  term  is  negligible. 


At  this  point  some  comments  about  the  values  of  the  functions  Tj_ (o^ ) 
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and  T2(o-j,o-j)  given  in  Figure  4.2  are  worthwhile.  Disregarding  the 
diagamma  function  (this  approximation  is  discussed  in  Chapter  5)  these 
functions  are  given  by: 


sinh(Bo.)  Bo.  Bo. 

W  =  cosh(Bo . )+cos(8a)  ~  tanh(“2~)  ~  ~ 

T  rn  a  i  -  sin(Ba) _ Ba  „  Ba 

’  v  cosh(Bo. )+cos(Ba)  cosh(Boi)+l  2 


(4.1?) 

(4.18) 


when  Ba  <<  1  and  Soi  «  1. 

The  values  of  T x ( a i )  and  12(0,, a,)  given  in  Figure  4.2  are  exactly 
reproduced  by  the  approximations  (4.17)  and  (4.18).  The  very  low 
value  of  Ti(o-j)  (=2.43*10*3)  indicates  that  at  this  low  resonance  wave 
number  stimulated  emission  is  an  important  process  with  regard  to  dis¬ 
persion  (and  absorption). 


Figure  4.3  illustrates  the  situation  when  the  resonance  wave  number  is 
1QQ  cm*1. 


Figure  4.3  Wave  number  dependence  of  the  different  terms  in  Equation 
(4.14)  for  resonance  wave  number  oi  =  100  cm*1 
The  functions  Tj(oi)  and  12(0,01)0=0.  have  the  values 
indicated  1 
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The  results  in  Figure  4.3  shows  that  the  peak  value  of  the  first  term 
in  Equation  (4.16)  is  dominating  with  at  least  two  orders  of  magnitude 
over  the  other  terms. 

The  values  of  T ^ ( o i )  and  T2(°i »ai )  stated  in  Figure  4.3  are  almost 
exactly  reproduced  by  the  approximations  (4.17)  and  (4.18).  The  value 
of  T i ( o -j )  (0.24)  has  increased  by  two  orders  of  magnitude  indicating 
the  reduced  importance  of  stimulated  emission  as  the  resonance  wave 
number  increases. 

The  part  of  the  approximations  in  (4.17)  and  (4.18)  which  are  valid 
when  only  0a  <<  1  shows  that  T j ( o i )  -*  1  and  T2(o,o-j)  -*  0  when  the  re 
sonance  wave  number  increases.  For  >  500  cm-1,  T i  ( o -j )  >  0.84  and 
7 2 ( 0 . o i )  <  8.5-10*5  Wften  the  temperature  is  296  K. 

As  a  result  of  this  discussion  it  is  reasonable  to  assume  that,  at 
standard  atmospheric  pressure  (Lorentz  broadning  dominates)  and 
temperature,  the  dispersion  in  the  infrared  and  visible  part  of  the 
spectrum  is  well  approximated  by 


An(o)  = 


1 

4tP 


l 

ij 


(mj) 


Si  (T) 


(v°)! 


-nr 


(4.19) 


The  result  in  Equation  (4.19)  could  easily  be  obtained  by  doing  the 
following  approximations  in  the  expression  for  the  optical  depth 
(Equation  (3.1))  tanh(0o/2)  =  1  and  f(-o,o^)  =  0  prior  to  the  Hilbert 
transformation. 


A  general  line  by  line  calculation  of  the  dispersion  valid  for  differ 
ent  pressures  and  temperatures  may  in  principle  be  implemented  by 
means  of  the  Equations  (4.5),  (4.7)  and  (4.13).  However,  the  presence 
of  the  digamma  function  in  Equation  (4.5)  makes  such  an  implementation 
rather  awkward.  In  the  next  chapter  approximations  that  simplify  the 
expression  of  Equations  (4.5)  and  (4.7)  will  be  made.  The  validity  of 
these  approximations  in  different  domains  of  wave  numbers  and  thermo¬ 
dynamic  parameters  will  be  discussed.  One  of  these  approximations 
will  be  implemented  as  a  computer  program  for  line  by  line  calculation. 
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5  APPROXIMATIONS  AND  IMPLEMENTATION  TN  FASC0D2 

The  aim  of  the  work  described  in  this  chapter  is  to  develop  an 
approximation  for  the  dispersion  formulas  calculated  in  Chapter  3 
whicn  is  suitable  for  line  by  line  calculation  within  the  framework  of 
FASCOD2  (10). 

The  first  part  of  this  chapter  will  mainly  discuss  the  importance  of 
the  different  terms  in  the  dispersion  formulas  especially  the  LorenLz 
dispersion  function  given  by  Equations  (4.9)  and  (4.5).  The  result  of 
this  discussion  will  establish  certain  domains  for  the  resonance  wave 
numbers  and  tne  the.'modynami c  parameters  where  different  approxima¬ 
tions  for  the  Lorentz  dispersion  function  are  valid. 

The  second  part  of  this  chapter  describes  one  specific  approximation 
and  the  implementation  of  this  approximation  in  FASCOD2  for  line  by 
Mne  calculation  of  dispersion. 


5- 1  Approxima1- ions  for  the  Lorentz  dispersion  function 

The  discussion  in  this  part  will  concentrate  on  approximations  for  the 
Lorentz  dispersion  function,  F2(o,Oi),  given  by  Equation  (4.9)  and 
(4-5).  The  reason  is  that  a  valid  approximation  *or  the  Lorentz  dis¬ 
persion  function  also  imply  a  valid  approximation  for  the  Voigt  dis¬ 
persion  function  given  by 


(0,0^ 


1  +7  x'°i  * 

—  /  F£(o,x)  •  exp[-ln2  (— -)  ]dx 

d  -®  d 


(5.1) 


This  statement  requires  that  the  exponential  in  Equation  (5.1)  has 
diminishing  values  near  the  edge  of  the  validity  reqion  for  the 
approximate  Lorentz  dispersion  function.  In  other  wo, ds  the  vali¬ 
dity  region  for  the  approximate  Voigt  dispersion  function  will  be 
some  Doppler  half  widths  less  than  the  validity  region  for  the 
approximate  Lorentz  dispersion  function. 
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from  Equation  (4.51,  (4.9)  ana  (4.14)  the  complete  expressions  for 
the  Lorentz  dispersion  functions  are 


i  °r° 


^j(°’°i)  -  ti  (o  -o)2+aJ  '  ^  1  ^ 0  i  ^  4  n  (a . -o) 2+oJ  ^(0,0. ) 


(5.2) 


O  40 


■  i  •  W  * ;  •  Tv°W  •  V3-0,’ 


where 


Tl(oi> 


2(0’V 


sinh(go-) 


rHhTFo’i'Ico^rai^i1"1  H<v1a>>-«,4+H<v,0») 


IB, 


lea) _ 

30-T+C0S 


(5.3) 


'  i 


^«v>a))+-K«4S(oj+,a)i! 


*  iRe  Ji°)l 


In  order  to  do  a  line  by  line  computation  of  dispersion  based  on  the 
Voigt  dispersion  function.  Equation  (5.1),  or  even  the  Lorentz  disper 
$iun  function,  the  functions  Tj(oi)  and  T2(o,o-j)  are  rather  awkward. 
Especially  the  digamma  functions  with  complex  arguments. 

It  is  therefore  necessary  to  analyse  the  value  of  the  terms  in  Eqja 
lion  (5.3)  in  order  to  determine  a  domain  of  resonance  wave  numbers 
and  t.crmodynami c  parameters  where  the  digamma  functions  can  be 
neglected.  This  analysis  is  described  in  detail  in  Appendix  B.  The 
result  is  given  in  Equation  (B.10)  stating  that  the  diagamma  func¬ 
tions  can  be  neglected  when 

sinh(3o  )  3 

cosh(0o.  )+cos(0a)  >>  if7 

Since  the  value  of  B*a  (=  hc/kT*a)  normally  is  much  less  than  one  (in 
atmospheric  physics  B*a  <  10“ * )  the  inequality  (5.4)  is  equivalent  to 


26 


o  >>  |  arctanh  (^y  6a)  =  §y  a  (5.5) 

in  atmospheric  physics  the  maximum  value  of  the  Lorentz  half  width  is 
of  the  order  0,5  cm-1.  A  neglection  of  the  diganma  functions  is  then 
a  valid  approximation  when 

oi  >>  0,3  cm'1  ■  (5.6) 

This  means  that  the  digamma  functions  have  some  influence  on  atmos¬ 
pheric  dispersion  for  resonance  transitions  in  the  wavenumber  interval 
up  to  -2  cm'1.  In  the  visible  and  infrared  wave  number  domain  the 
influence  of  these  functions  are  of  no  importance,  but  atmospheric 
influence  on  electromagnetic  radiation  with  wave  number  less  than 
2  cm'1  ( correspondi ng  to  a  frequency  of  30  GHz)  is  of  great  importance 
in  many  technological  applications. 

When  the  inequality  (5.4)  holds  it  is  possible  to  define  approximate 
Voigt  and  Lorentz  dispersion  functions,  frvA(°«°i)  and  fr4A(°»°i) 

f  v  A ( 0 ' 0 1  >  =  /  'raA(o.x)-fd(x,o1)dx  (5.7) 

where 

fd(x,o.)  -  The  Doppler  line  shape  function 


1  v° 


1  a 


2A^0,0i  ^  "  it  (o^  -o)  J+aJ  ^  1 A  ^ 0  i  ^  +  i  (o  ,  -o)  s+aJ  T2A^°i^ 


(5.8) 


TlA<°i> 


sinh(8o.) 


cosh(6o1 )+cos(8a) 


(5.9) 


T2A<°i> 


sin(Ba) 


coshCBo^+cosCBa) 


In  order  to  get  a  picture  of  the  importance  of  the  different  terms  and 
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variables  in  the  approximate  Voigt  dispersion  function  a  computer 
program  that  performs  the  integration  in  Equation  (5.7)  for  one 
single  resonance  transition  is  developed.  In  this  case  the  dispersion 
is  given  by  the  simplified  formula 

An(o)  =  ^  W(m.)  S.  ^(o.o^  +  FvA(o,-0i)]  (5. If) 


The  programme  calculates  only  the  sum  of  the  two  approximate  Voigt 
dispersion  functions  in  Equation  (5.10).  This  sum  can  be  expressed  as 
a  sum  of  a  Dominating  term,  FvAD(°>°i)»  anci  a  Correcting  term, 
FvAC(°-°i).  v.i.z 


Fva(°-01>  +  FvA<C--ai>  '  FvAD<°'°i>  +  FvAC(0'°1> 


(5.n; 


where 


FvAD(0’°i)  =  \  T1A(x)  I  (x-a)*+a2  Vx,0i)dx 


(5.12) 


FvAC(0^°i)  =  [  T1A(x)  i  (x+o) *+a2  fd(x,oi)dx 


+.i  r2A(x)  *  *  (x-o°2+o2 


^  (x+o)2+a2  *  fd(x,0i)  dx 

(5,13) 


f  (x  0  )  =  J  .  — 

d'*’  1 '  '  it  a . 


x-0i 

exp[-ln2  (-—l)2] 


(5.14) 


The  calculation  of  the  sum  in  Equation  (5.11)  and  the  correcting  term 
in  Equation  (5.13)  is  done  for  two  different  thermodynamic  situation. 
One  in  the  high  pressure  limit  where  the  homogeneous  (Lorentz)  broad¬ 
ening  dominates  and  one  in  the  Intermediate  regime  where  the  homogene¬ 
ous  and  Inhomogeneous  (Doppler)  line  widths  are  equal.  The  calcula¬ 
tions  are  done  for  two  resonance  wave  numbers  *  1  cm-1  and 
Oi  =  100  cm-1.  The  Doppler  line  width  used  in  the  calculations  is 
adequate  for  a  H£0  molecule. 
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The  Figures  5.1  and  5.2  show  both  the  sum  of  the  dominating  and  cor¬ 
recting  term  and  the  correcting  term  only  for  the  two  different  ther¬ 
modynamic  situations  when  the  resonance  wave  number  is  1  cm-1. 


Figure  5.1  The  total  Voigt  dispersion  function  (5.11)  and  the  cor¬ 
rect  i rig" Term'  (5.12)  at  standard  temperatureand  pressure 
when  Oi  s'  1  cm"1 

(6  =  4.86- IQ"5  cm,  Tia(oi)  =  2.43-10"3,  T2a(oi)  =  2.43*10“M 


Figure  5.2  The  total  Voigt  dispersion  function  (5.11)  and  the  cor¬ 
recting  term"  (5.12)  when  the  Lorentz  and  Doppler  half 
widths  are  equal  and  oi  *  1  cm-1 

(0  «  6.23- MT’ "cm,  T 1 A ( 0 1 )  «  3,11*  10"3  ,T2a(°1  )  a  3-96-10"  •’J 
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The  big  differences  between  the  shapes  and  peak  values  in  Figure  5.1 
and  5.2  are  mainly  due  to  the  five  order  of  magnitude  change  of  the 
effective  line  width.  This  is  also  the  reason  for  the  vanishing  cor¬ 
recting  term  in  Figure  5.2.  When  the  results  in  Figure  5.1  and  5.2  are 
used  for  calculation  of  the  dispersion  by  Equation  (5.10)  (assuming 
equal  line  strength  S^)  the  number  density  must  be  scaled  properly  by 
the  factor  (p/Po) * (T0/T) ) .  The  resulting  peak  dispersion  is  almost 

equal  for  the  two  situations. 

Figure  5.3  and  5.4  show  almost  the  same  situations  as  the  preceeding 
figures  when  the  resonance  wave  number  is  100  cm-1. 

ine  difference  in  shape  and  peak  values  between  Figure  5.3  and  5.4  is 
mainly  due  to  the  tree  order  of  magnitude  change  of  the  effective  line 
width.  The  correcting  term  vanishes  in  both  situations  mainly  because 
the  resonance  wave  number  is  increased  by  two  orders  of  magnitude. 

The  peak  value  in  Figure  5.4  scaled  by  the  number  density  factor  is 
almost  equal  to  the  peak  value  in  Figure  5.3  resulting  in  about  equal 
peak  dispersion  for  the  two  situations  (assuming  equal  line  strength 

Si). 

The  peak  value  of  the  dispersion  for  resonance  wave  number  100  cm-1  is 
however  about  two  orders  of  magnitude  greater  than  the  peak  value  of 
the  dispersion  for  resonance  wave  number  1  cm-1.  This  is  seen  by  com¬ 
paring  the  peak  values  In  Figure  5.1  and  5.3.  (The  same  result  is 
obtained  by  comparing  the  peak  values  of  Figure  5.2  and  5.4  after 
proper  scaleing  by  the  number  density  factor).  The  reason  is  that  the 
value  of  the  function  Tia(oi)  is  almost  linear  with  wave  number  in  this 
region.  This  again  demonstrates  the  reduced  importance  of  stimulated 
emission  with  increasing  resonance  wave  number. 
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Figure  5.3  The  total  Voigt  dispersion  function  (5.11)  and  the  cor¬ 
recting  term  (5.12)  at  standard  temperature  and  pressure 
when  oi  =  100  cm-1 

(0  =  4.86- lO"3  cm,  Tia(oi)  =  0.24,  T2a(°i)  ■  2.31-10"*) 


Figure  5.4  The  total  Voigt  dispersion  function  (5.11)  and  the  cor 
rectlnq  term  ;5. 12)  when  the  Lorentz  and  Doppler  haTf~ 
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5.2  Implementation  in  FASC0D2 

The  discussion  in  this  part  will  concentrate  on  the  implementation  of 
an  appropriate  formula  for  line  by  line  computation  of  dispersion  by  a 
sliahtlv  modified  FASC0D2. 


For  line  by  line  computation  of  absorbtion  the  Voigt  line  shape  func¬ 
tion  of  Equation  (3.5)  is  used.  This  necessitates  a  numerical  inte¬ 
gration  to  achieve  each  value  of  the  Voigt  line  shape.  To  speed  up 
this  integration  process  a  special  algorithm  has  been  developed  based 
on  the  fact  that  the  Voigt  line  shape  has  a  "fall  off"  -  (o-o^)‘2  div 
tant  from  the  line  center  while  the  inhomogenous  (Doppler)  broadning 
is  important  near  the  line  center. 


This  algorithm  will  probably  not  give  a  correct  result  for  the 
approximative  Voigt  dispersion  function  of  Equation  (5.7)  since  the 
dominating  term  has  a  "fall  off"  ~  (o-Oi)'1  (see  Equation  (5.12))  dis 
tant  from  the  line  center.  To  implement  the  Voigt  dispersion  function 
in  FASC0D2  a  new  algorithm  for  numerical  integration  has  to  be 
developed.  This  problem  is  not  addressed  in  this  report. 


Due  to  the  lack  of  a  new  numerical  integration  algortihm  only  the 
approximative  Lorentz  dispersion  function  has  been  implemented  in 
FASC0D2.  The  implemented  formula  for  line  by  line  computation  of 
dispersion  is 


An(°)  »  l  WOnj)  $i (T)  l>:iA(°.°i)+FaA(o,-oi)] 
where 


o,~o 


1  i  1  a 

F£A^0,0i^  "  W  (o^-o^+o2  TlA(oi)+i  (o^-o)2+a2  T2A^°i^ 

1  °i+0  1  a 

P£A^0’”°1^  “  i  (o.j  +o ) 2  +o2  T1A  ^  ^  +tt  (o^+o)  2+a2  T2A^°i^ 

s1nh(goi ) 

TlA<°i>  =  cosh(0o.j )+cos(0a) 


(5.15) 


(5.16) 


(5.17) 


(5.18) 
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t  (a  )  -  _ 5^n(^g) _  (5.19) 

2A'°i'  cosh(So^  )+cos(0a)  '  ' 

3  =  hc/kT 

a  -  Lorentz  half  width  (HWHM) ;  function  of  temperature  and  pressure 

The  connection  between  the  line  strength,  Si (T) ,  in  Equation  (5.15)  and 
the  line  strength,  Si ( T) ,  filed  on  the  HITRAN  molecular  database  (11)  i 

S  (T) 

vt>  - — Ho”  (5‘2°) 

oitanh(~) 

The  units  of  the  HITRAN  line  strength  is  cm/molecule. 

Equation  (5.15)  should  give  an  adequate  description  of  the  dispersion 
in  a  mixture  of  gases  when  the  following  conditions  are  fulfilled: 

a)  The  dipole  approximation  and  a  classical  description  of  the 
radiation  field  are  valid 

b)  Line  coupling  effects  are  negligible 

c)  Homogenous  (Lorentz)  broadning  is  dominant 

d)  The  resonance  wave  number  0}  >  1  cm-1 

The  theoretical  description  in  this  report  makes  it  possible  to  defeat 
condition  d)  by  introducing  the  digamma  functions  in  Equations  (5.18) 
and  (5.19),  and  condition  c)  by  developing  a  new  numerical  integration 
algorithm  for  the  Voigt  dispersion  function. 
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6  LINE  BY  LINE  CALCULATION  OF  ATMOSPHERIC  DISPERSION 


This  chapter  describes  results  from  some  line  by  line  calculations  of 
atmospheric  dispersion  by  FASCOD2  modified  according  to  Equations 
(5.13;-(5.19).  Four  results  are  presented  and  discussed.  These 
results  show  the  atmospheric  dispersion  in  different  wave  number 
regions  at  different  spectral  resolution  for  the  Midlatitude  summer 
model,  i  e  temperature  294  <,  pressure  1013  mb  and  a  water  content  of 

3 

14  g/m  at  ground  level . 

Figure  6.1  shows  a  result  where  the  contribution  to  the  dispersion 
steems  mostly  (though  not  totally),  from  one  single  transition  in  the 
water  molecule  with  resonance  wavenumber  oi  =  1684.837  cm-1  and 
Lorentz  half  width  a  =  0.0962  cm-1. 
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Figure  6.1  The  atmospheric  dispersion  at  ground  level  in  the  wave 
number  Interval  1684-1686  cm~l  for  the  Midlatitude  sum¬ 
mer  model. 

The  main  contribution  to  the  dispersion  is  from  a  transi¬ 
tion  in  the  water  molecule  with  resonance  wavenumber 
Oj  *  1684.837  cm-1. 


The  result  In  Figure  6.1  shows  the  classical  dispersion  connected  to 
one  single  resonance  transition.  The  maximum  and  minimum  value  of  the 
dispersion  in  Figure  6.1  are 
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An  =  3.9-10-*  (6.1) 

max 

An  =  -1.0- 10-*  (6.2) 

mi  n 

For  this  single  line  in  the  infrared  part  of  the  spectrum  it  is  easy 
to  verify  the  result  in  Figure  6.1  from  Equation  (4.19).  For  a  single 
line  Equation  (4.19)  is  modified  to  read 


An(o) 


W(m1 ) ♦  Si(T)- 


To 


0 

+  a5 


(6.3) 


By  introducing  the  HITRAN  line  strength  from  Equation  (5.20)  and  the 
approximation  tanh(0o. /2)=1  in  the  infrared  at  standard  temperature. 
Equation  (6.3)  is  modified  to 


o4-  o 


An(o 


>  •  JF  •  5;  •  ;o1-’o)’  +  a;  •  “<"V  •  VT> 


(6.4) 


The  maximum  and  minimum  value  of  Equation  (6.4)  occur  when  the  wave 
number  is  equal  to  oi  ±  a  and  the  peak  values  of  the  dispersion  are 


.  1  W("V  •  VT> 

An  =  ±  = — r  •  - 

max, min  Sir  o.  •  a 

when 


(6.5) 


o 


max, min 


a 


(6.6) 


The  number  density  of  water  molecules  and  the  line  strength  on  the 
HITRAN  data  base  for  this  transition  gives  the  following  value  for  the 
product  of  number  density  and  line  strength 


W(mi)  •  St(T)  =  0.1304  cm"2  (6.7) 

Wave  number  and  value  for  the  peak  dispersion  calculated  from  Equation 
(6.5)  and  (6.6)  is 
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Anmav  =  ±  1.0* 10~s 
max ,mi n 


o  ,  =  1684.741  cm*1 
max 


om.„  =  1684.933  cm"1 
mi  n 


(6.8) 


(6.9) 


The  calculated  result  in  (6.8)  and  (6.9)  are  in  very  good  agreement 
with  the  results  presented  in  Figure  6.1.  Wave  number  and  value  for 
the  minimum  dispersion  have  the  best  agreement  due  to  less  influence 
from  another  resonance  line  at  lower  wavenumbers  in  the  FASC0D2 
calculations. 


The  conclusion  from  the  preceeding  comparison  between  the  FASC0D2 
result  and  a  theoretical  result  is  that  the  implementation  in  FASC0D2 
most  probably  performs  a  correct  calculation  of  the  dispersion. 

Figure  6.2  shows  the  atmospheric  dispersion  in  the  wavenumber  region 
35-45  cm'1,  where  several  transitions  in  the  water  molecule  are  con¬ 
tributing  to  a  rather  complex  result  due  to  interference  between  dif¬ 
ferent  lines.  The  atmospheric  model  is  still  Midlatitude  summer  at 
ground  1 evel . 
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Figure  6.2  The  atmospheric  dispersion  at  ground  level  in  the  wave 
number  interval  35-45  cm"1  for  the  Midlatitude  summer 
model . 

Several  transitions  In  the  water  molecule  with  different 
line  strength  are  contributing. 
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The  complex  nature  of  the  dispersion  in  Figure  6.2  is  due  to  the 
interference  between  a  number  of  different  water  vapour  lines  with 
different  line  strength.  The  peak  dispersion  in  Figure  6.2  is  about 
one  order  of  magnitude  greater  than  the  peak  dispersion  in  Figure  6.1. 
This  effect  can  be  infered  from  Equation  (6.5)  although  this  equation 
is  not  strictly  correct  in  this  wavenumber  region. 

Figure  6.3  shows  the  atmospheric  dispersion  for  the  entire  wavenumber 
region  0-500  cm-1.  The  majority  of  the  "lines"  are  due  to  transitions 
in  the  water  molecule. 
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Figure  6.3  Atmospheric  dispersion  at  ground  level  in  the  wave- 
number  interval  0-500  cm" 1  for  the  Midlatitude  summer 
model . 

The  main  contribution  to  the  dispersion  steems  from  seve¬ 
ral  transition  in  the  water  molecule 


The  maximum  and  minimum  value  of  the  dispersion  in  Figure  6.3  is 
1.06* 10-3  for  the  resonance  line  near  150  cm~l.  By  anticipating  the 
validity  of  Edlen's  formula,  see  (1),  for  dry  air  in  the  p-wave  region 
it  is  possible  to  calculate  the  maximum  and  minimum  value  of  the 
refractive  Index  around  ~150  cm-1  for  the  gas  mixture  in  this  example. 
Edlen's  formula  corrected  for  a  temperature  of  294  K  gives  (o  =  150  cm'1)* 
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(n-l)dpy=  2.6*10"*  (6.10) 

The  influence  of  water  vapour  on  the  refractive  index  then  gives  the 
following  maximum  and  minimum  value  for  the  refractive  index  of  the 
gas  mixture  near  the  wave  number  150  cm"1 


max 


nmi  n 


1  +  0.26-10" ’  ±  1.06-10"3 


'1.00132 

0.99920 


(6.11) 


The  calculations  in  (6.11)  -show  the  rather  peculiar  result  that  the 
minimum  refractive  index  is  less  than  unity.  The  calculations  o 
show  that  every  resonance  line  in  Figure  6.3  with  a  minimum  'a  1 u  for 
the  dispersion  less  than  2.6*10"*  will  exhibit  a  corresponding  minimum 
refractive  index  less  than  unity.  The  majority  of  the  resonance  lines 
in  this  wave  number  region  will,  as  seen  from  Figure  6.3,  exhibit  this 
pecuiarity.  Whether  these  results  are  correct  remains  to  be  deter¬ 
mined  by  experiments. 


Figure  6.4  shows  the  atmospheric  dispersion  in  the  wave  number 
interval  2300-2400  cm"1  for  the  Midlatitude  summer  model,  at  ground 
level.  It  is  mainly  the  CO2  molecule  that  contributes  to  the  disper¬ 
sion  in  this  wavenumber  interval. 


The  structures  of  the  dispersion  in  Figure  6.4  are  due  to  the  P-  and 
R-branch  of  viberational -rotational  transition  of  CO.,  centered  at  wave 
number  2349.146  cm"1. 

In  this  wave  number  interval  the  Edlen  formula  gives  the  following 
value  for  the  refractive  index  of  dry  air  (0  =  2350  cm"1) 

(n-l)dry  =  2.7*10"*  (6.12) 

The  peak  dispersion  in  Figure  6.4  are  about  two  orders  of  magnitude 
less  than  the  Edlen  value.  Resonance  lines  in  this  wave  number  region 
therefore  contribute  negligible  to  the  refractive  index  of  air  in  the 
Midlatitude  summer  model. 
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Figure  6*4  Atmospheric  dispersion  at  ground  level  in  the  wave  number 
interval  2300-2400  cm-1  for  the  Midlatitude  summer  "model  7 
It  is  matr.ly  CO?  that  contributes  to  the  disperison  in 
this  interval. 


An  overview  of  the  dispersion  for  the  entire  wave  number  region 
500  cm-1  -  15000  cm-1  for  the  Midlatitude  summer  model  is  given  in 
Appendix  C. 

The  input  data  for  FASC0D2  are  as  mentioned  provided  by  the  line  data 
base  HITRAN.  This  database  contains  l'ne  data  for  molecular  transi¬ 
tions  with  resonance  wave  numbers  in  the  wave  number  interval 
0-17900  cm-1.  In  the  wave  number  interval  aoove  5000  cm'1  the  disper¬ 
sion  is  a  very  small  quantity  in  air  (less  than  3*10-7  for  the  Mid¬ 
latitude  summer  model).  It  has  not  been  possible  to  find  measurements 
with  an  acruracy  that  makes  the  verification  of  the  fine  structure  of 
our  computations  possible  in  this  Interval.  Below  5000  cm-1  there  are 
wave  number  1ntervals  where  the  dispersion  has  a  substantially  higher 
value.  The  most  interesting  domain,  when  it  comes  to  comparl.Rns 
between  calculation  and  measurements,  Is  the  wave  number  interval 
0-500  cm-1.  In  this  wave  number  Interval  the  modified  version  of 
FASC002  predicts  that  the  value  of  the  disperison  Is  In  the  range 
±1.2* 10-3  for  the  Midlatitude  summer  model. 


DiSPETiiON  2?  DISPERSION 
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A  paper  by  Liebe,  (3),  investigates  a  model  made  for  the  computation 
("c  attenuation  and  dispersive  delay  in  the  wave  number  interval 
j-33  cm-1.  Comparisons  between  his  model  and  measurements  are  done* 
anc  Liebe  ciaims  that  good  agreement  between  model  predictions  and 
field/laboratory  experiments  is  obtained. 

Figure  5.5  shows  trie  results  from  Liebe' s  model  for  the  dispersion, 
whilst  Figure  6.6  shows  the  dispersion  calculated  by  the  modified  ver¬ 
sion  of  FASC0D2.  The  gas  mixture  and  thermodynamic  parameters  are 
identical . 


Qure  6.5  Ihe  calcualted  dispersion  for  air  where  the  relative 
humidity  is  lOOTi.  the  temperature  is  303  <  and  the 
pressure  is  1013  mb.  From  Liebe.  (3)7 
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Figure  6,6  The  dispersion  calculated  by  the  modified  version  of 

FASC0D2.  - 

The  gas  mixture  and  the  thermodynamic  parameters  are 
identical  to  those  used  by  Liebe  (3). 
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A  close  comparison  between  the  dispersion  in  Figure  6.5  and  6.6  shows 
excellent  agreement  especially  between  the  peak  dispersion  at 
18.5  cm-1  wave  number. 

The  results  presented  in  this  chapter  demonstrate  that  the  modified 
version  of  FASC0D2  performs  a  correct  line  by  line  calculation  of  the 
dispersion  due  to  resonance  transitions  when  homogenous  (Lorentz) 
broadening  is  dominant. 


7  CONCLUSIONS  AND  PROPOSITIONS  FOR  FURTHER  RESEARCH 

From  the  connection  between  the  complex  propagation  constant  of  an 
electromagnetic  field  in  a  medium  and  the  complex  linear  suscepti¬ 
bility,  a  Hilbert  transform  relationship  between  the  anomalous  dis¬ 
persion  and  the  optical  depth  has  been  established  via  the  Kramer- 
Kronig  relations.  Based  on  the  Van  Vleck-Huber  formalism  for  the 
optical  depth  an  explicit  expression  for  the  dispersion  in  a  mixture 
of  gases  is  developed.  This  expression  determines  the  dispersion 
whether  the  line  shape  is  described  by  the  Voigt,  Lorentz  or  Doppler 
1 ine  shape  functions. 

The  derived  dispersion  is  valid  for  all  wavenumbers  provided  the 
dipole  approximation  is  valid  and  a  classical  description  of  the 
electromagnetic  field  is  adequate.  The  effect  of  line  coupling  is  not 
taken  into  account.  The  Hilbert  transformation  of  the  optical  depth 
is  performed  without  approximations. 

The  importance  of  the  different  terms  in  the  derived  expression  is 
visualized  and  discussed  by  numerical  examples. 

Since  the  complete  expression  for  the  dispersion  is  rather  awkward  for 
numerical  computations  an  approximation  is  developed.  This  approximate 
expression  for  the  dispersion  Is  valid  when  the  resonance  wave  number 
for  a  transition  is  greater  than  1  cm-1  (30  GHz).  Numerical  examples 
of  the  approximate  dispersion  for  a  single  resonance  line  described  by 
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a  Voigt  line  shape  function  is  also  visualized  and  discussed. 

In  order  to  perform  a  line  by  line  computation  of  the  dispersion  in  a 
mixture  of  gases  the  computer  code  FASC0D2,  developed  by  the  US  Air 
Force  Geophysics  Laboratory  (AFGL) ,  has  been  modified.  The  modified 
version  of  FASC0D2  only  performs  correct  computation  of  the  dispersion 
when  the  gas  pressure  is  in  the  regime  where  homogenous  (Lorentz) 
broadning  is  dominant  and  the  resonance  wave  number  is  greater  than 
1  cm"1. 

Examples  of  line  by  line  calculated  atmospheric  dispersion  for  the 
Midlatitude  Summer  model  at  standard  atmospheric  pressure  are  given. 
One  example  shows  that  the  calculated  dispersion  from  FASC0D2  is  iden¬ 
tical  to  the  theoretical  result  for  a  specific  transition  in  the 
infrared  region.  This  result  confirms  that  the  modifications  in 
FASC0D2  are  correct.  Another  example  from  the  wave  number  region 
1-33  cm"1*  is  a  comparison  between  the  calculated  result  from  FASC0D2 
and  the  measurements  of  H  J  Hebe  at  the  Institute  for  Telecommunica¬ 
tion  Sciences,  USA.  These  two  results  are  in  excellent  agreement. 

The  result  from  these  two  examples  leads  to  the  conclusion  that  the 
theoretical  expression  for  the  dispersion  derived  in  this  report  and 
the  implementation  in  FASC0D2  gives  a  correct  line  by  line  computation 
of  the  dispersion  in  a  mixture  of  gases  when  line  coupling  effects  are 
negl igible. 

In  order  to  achieve  increased  confidence  in  the  theoretical  results 
derived  in  this  report  a  greater  effort  is  needed.  The  rest  of  this 
chapter  is  therefore  devoted  to  propositions  for  further  research. 

Results  from  the  existing  modified  FASC0D2  should  be  compared  with 
high  resolution  measurement  of  dispersion  In  different  wavenumber 
regions  at  different  temperature  and  humidity. 

In  order  to  do  line  by  line  calculations  for  resonance  transitions 
<  1  cm-1  (30  GHz)  the  neglected  digamma  functions  in  the  lorentz  dis¬ 
persion  function  should  be  Incorporated  in  the  existing  line  by  line 
code.  This  would  facilitate  comparison  between  calculated  and  mea- 


42 


sured  dispersion  for  frequencies  below  30  GHz;  a  frequency  range  of 
great  technological  importance. 

To  get  a  complete  description  of  the  dispersion  valid  for  all  gas 
pressures  the  Voigt  dispersion  function  should  be  used  by  the  computer 
code.  This  requires  the  development  of  a  fast  numerical  integration 
algorithm  for  the  calculation  of  the  Voigt  dispersion  function. 

The  hypothesis  that  the  refractive  index  of  air  in  the  visible  and 
infrared  is  the  result  of  adding  the  "distant  wing  dispersion"  due  to 
ail  strong  UV-transi tions ,  should  be  tested  by  the  existing  modified 
FASC0D2.  This  requires  that  data  for  the  relevant  UV-transitions 
exists  on  the  HITRAN  data  base,  and  a  change  in  FASC0D2  of  the  maximum 
distance  from  the  line  center  where  the  effect  from  a  line  is  taken 
into  account.  The  result  of  this  line  by  line  calculation  depends 
heavily  on  the  correctness  of  the  Lorentzian  "wings"  but  even  a  result 
of  the  same  order  of  magnitude  as  given  by  the  Edlen  formula  would  be 
an  encouragement. 
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APPENDIX  A 

DETAILS  OF  THE  HILBERT  TRANSFORMATION  IN  CHAPTER  4 


i  lm  {o'} 

0  (2n  +  1)n 

!  P 


Figure  A.l  The  integration  contour  C  in  the  complex  plane  which  is 
used  in  the  evaluation  of  the  integral  in  equation  (A.l) 
The  poles  of  the  integrand  of  equation  (A.l)  are  ' 

indicated  by  the  symbol  O 


In  Chapter  4,  Equation  (4.4),  it  was  shown  that  in  the  limit  where 
R  -  and  e  -*  0 


PP  J  F ( a 1 ) do 1 


2th 


fThe  sum  of  the  residues  of 
* F ( o ■ )  inside  the  contour  C 


}  -  1 im  /  F(o' )do 
e-*o  P 
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where 


F(o') 


tanh  (|2j 


1  1  a 

o' -o  u  (o'-x)z+gJ 


*  tanh  (|^) 


1  1  a 

o'-o  i  o' -x-1a 


1 

o' -x+ia 


(A.2) 


The  integration  of  the  function  in  Equation  (A.2)  along  the  curve  f 
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can  be  evaluated  as  follows.  The  pole  of  F(o')  at  o'=o  is  simple,  so 
that  the  Laurent  series  has  the  form 


M°)  ;  n 

F(o')  =  fa  + 


By  doing  the  substitution 
o'-o  =  e  ei0 

the  integral  along  P  is  given  by 


0  b.(o)  iS 
/  F ( o 1  )do*  =  I  1  i  — 

P  it  ee 


.18 


ieel8d0  +  X  a  .  en  i  /  ein  d8 
n=l  tt 


*  -”1  »1<°>  c"  ^n"3  3 


In  the  1  imit  e-»0  one  gets 

lim  /  F(o' )do‘  =  -ni b, (o) 
e->o  I” 


(A. 3) 


(A. 4) 


(A. 5) 


(A. 6) 


bi(o)  is  the  residue  of  F(o')  at  o'=o.  From  Equation  (A. 2)  this 
residue  is  given  by 


bj(o)  =  res  F(o')  =  tanh  (|^)  i  (q_x)»+q»  {A*7) 

In  the  limit  e-»o  the  Integral  along  P  is  then 

lim  /  F(o'  )do‘  =  -ir1*tanh(|2)  •  i  •  (A-8) 

To  evaluate  the  integral  in  Equation  ( A. 1 )  the  sum  of  the  residues  of 
F(o')  inside  contour  C  remains  to  be  determined. 


One  of  the  poles  inside  the  curve  C  Is  a  simple  pole  at  o'  =  x+ia. 
The  residue  of  F(o')  at  this  point  is 


res  F(o')  =  tanh 
o'=x+1a 


-i  1  1 

J~  tt  (x-o)+io 


(A. 9) 


The  rest  of  the  poles  of  the  function  F(o')  in  Equation  (A. 2)  inside 
the  closed  contour  C  lie  on  the  imaginary  axes.  These  poles  are 
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infinitly  many,  simple,  and  they  steem  from  the  factor  tanh  (Bo')/2  in 
F(o’).  This  factor  is  singular  when 


^  =  i  \  (2n+l) , 


n  =  0,  1,  2, 


The  residue  of  pole  number  n  is 


res  F(o ' )  =  ~ 


1 


.  j2n+l 

o'  =  — n 


^  i  ~g—~TT-o  ( i~^~~^T~x) 2 +aJ 


(A. 10) 


(A. 11) 


Note  that  this  sequence  of  residues  goes  as  n~3  for  big  n,  so  that  the 
corresponding  sum  of  residues  converges.  To  find  an  explicit  expres¬ 
sion  for  the  series  one  can  use  partial  fraction  decomposition  and  the 
relation  (12) 


j0  (n+a)  (n+b)  ‘ 

where  ^(z)  is  the  digamma  function  which  is  defined  by 

V(z)  =  ^  (in  T(z))  =  (A.  13) 

r(z)  is  the  gamma  function  (9). 

The  sum  of  the  residues  on  the  imaginary  axis  is  after  some  algebra 


I  res  F(o')  = 
n=o  0.=i2ntin 


1 

Sr 


1  l 
n  (o-x)+ia 


*  if  (x-ia)) 


1 

IT 


+ 


1 

IT 


2ia 

(o-x)  2 -K3 J 


♦  (#  +  ^  °) 


(A. 14) 


By  using  the  Equations  (A. 8),  (A. 9)  and  (A. 14)  the  result  of  the 
Integration  Is 
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(A. 15) 


The  principal  part  of  the  integral  evaluated  is  by  definition  a  real 
function,  but  a  brief  look  at  the  right  hand  side  of  (A. 15)  can  hardly 
convince  anyone  that  this  is  so.  By  rewriting  (A. 15)  using  the  fol¬ 
lowing,  relations  (13) 

Imj4i(Vi+iy) }  *  tor  tanh  (Try)  (A. 16) 

4/(l-z)  -  'V(z)  =  tt  cotg  ttz  (A. 17) 

vjj(Z*)  =  i|j*(z)  (A. 18) 

the  result  in  (A. 19)  is  obtained  after  some  cumbersome  algebra 
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(A. 19) 


The  result  In  (A. 19)  Is  obviously  a  real  function. 
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APPENDIX  B 


NEGLECTING  THE  DIGAMMA  FUNCTIONS 


In  order  to  reduce  the  computational  burden  in  the  computation  of  the 
Lorentz  dispersion  function  it  is  necessary  to  establish  an  approxima¬ 
tion  where  the  digamma  functions  can  be  neglected.  The  Lorentz  dis¬ 
persion  function  defined  by  Equation  (5.2)  is: 


i  V  0 

it  o.-  o)s+a5 


sinn(8o. ) 

cos'hTsg^T+cos(3g)  *  i  ■"'[t(^(or>a))-«>4 


+  -  a 
+  u  (o.-  o)J+aJ 


-  |  Re{<p(%  +  g  a)} 


(BA) 


The  desired  approximation  is: 

1  c^-  o  sinh(3o. ) 

F2A^0,0i^  “  i  (o^-  a)1  +  a*  cosh(Boi )  +  cos(Ba) 

.  1  a  sin(ga) _  ,R 

tt  ( o . -  o)z  +  a2  cosh (30^ )  +  cos(Ba)  '  *  ' 

Note  that  the  peak  value  of  the  second  term  in  equation  (B.2)  will 
only  be  of  the  same  order  of  magnitude  as  the  peak  value  of  the  first 
term  when  the  resonance  wave  number  is  small.  However  the  wave  number 
dependence  of  the  second  term  is  such  that  it  will  always  be  unequal 
to  zero  when  the  first  term  equals  zero.  This  is  why  the  second  term 
is  kept  for  all  resonance  wave  numbers. 

The  purpose  of  this  appendix  Is  to  find  the  domain  of  resonance  wave 
number  and  thermodynamic  parameters  where  the  approximation  (B.2)  is 
valid.  In  order  to  compare  Equation  (B.l)  with  Equation  (B.2)  the 
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factors  containing  the  digamma  functions  in  (B.l)  are  expanded  in 
Taylor  series. 
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The  Figures  (B.l)  and  (B.2)  show  that 
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This  means  that 
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for  all  resonance  wave  numbers  o^. 

The  second  term  in  the  first  [•]  in  Equation  (B.l)  is  then  negligible 
compared  to  the  first  term  when 
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0  s  ^  Re  {^(x  +  i  y)}x=!^  <  3  when  y  >  0  (B.8) 

which  implies  that 
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The  fourth  term  in  Equation  (B.l)  is  then  neglible  compared  to  the 
first  term  in  (B.l)  when 


o^-  a  sinh(8o.)  a 
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A  simplification  of  Equation  (B.10)  yields 
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sinh  (Bo.) 

cosh  (Bc^)  +  cos(Ba) 


(B.11) 


Since  the  Equations  (B.7)  and  (B.ll)  are  identical  the  digamma  func¬ 
tions  in  (B.l)  can  be  neglected  when  Equation  (B.ll)  is  satisfied. 
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APPENDIX  C 

AN  OVERVIEW  OF  ATMOSPHERIC  DISPERSION  IN  THE  WAVE  NUMBER  INTERVAL 
500-15000  cm-1 

This  Appendix  gives  an  overview  of  calculated  atmospheric  dispersion 
in  the  wave  number  interval  500-15000  cm-1.  The  atmospheric  model  used 
in  these  calculations  is  the  Midlatitude  summer  model,  (pressure 
1013  mb,  temperature  294  K). 

The  interval  (0-500)  cm-1  is  not  presented  here  because  Figure  6.1  cov 
ers  this  spectral  range.  The  interval  15000-17900  cm-1  is  not  shown 
because  the  dispersion  is  /ery  small  in  this  region  (<<10~*). 

Each  *igure  cove's  a  broad  interval  of  wave  numbers.  These  figures 
therefore  gives  only  an  impression  of  the  magnitude  of  the  dispersion 
in  the  different  wave  number  intervals.  The  dispersion  varies  abou+ 
four  orders  of  magnitude  throughout  the  wave  number  interval;  from 
about  l3's-10-\ 
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Figure  C.3  The  atmospheric  dispersion  for  the  wave  number  interval 
6000-9000  cm-1.  Midlatitude  summer 
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C.4  The  atmospheric  dispersion  for  the  wave  number  interval 
9000-12000  cm'1.  Midlatitude  summer 
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C.5  The  atmospheric  dispersion  for  the  wave  number  interval 
12tid)Q- 15000  cm-1.  Midlatitude  summer 


